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Abstract 

We propose a new method to obtain the K vrvr amphtude from K ^ tt which 
allows one to fully account for the effects of final state interactions. The method is 
based on a set of dispersion relations for the K — > vrvr amplitude in which the weak 
Hamiltonian carries momentum. The soft pion theorem, which relates this amplitude 
to the K — > TT amplitude, can be used to determine one of the two subtraction 
constants - the second constant is at present known only to leading order in chiral 
perturbation theory. We solve the dispersion relations numerically and express the 
result in terms of the unknown higher order corrections to this subtraction constant. 



1. The K TTTT amplitude has been the subject of numerous studies over the years, the 
main challenge being either to explain the A J = 1/2 rule or to provide a reliable estimate 
of e' / e 0, 3 . Among these attempts, the lattice approach is in principle the most rigorous 
as the weak matrix elements are calculated from first principles in a truly nonperturbative 
way. However, the inclusion of final state interactions is problematic. The calculation of 
the K ^ TTTT amplitude, e.g., proceeds by calculating the K ^ n amplitude on the lattice 
and using chiral perturbation theory (CHPT) at tree level to obtain the physical decay 
amplitude [Q. As is well known, the latter step induces a sizable uncertainty in the final 
result, commonly estimated to be around 30%, the typical size of next-to-leading order 
corrections in chiral SU(3). To discuss the relation between the K ^ mr and the K it 
amplitude, it is necessary to allow the weak Hamiltonian to carry momentum. Then the 
former amplitude becomes a function of the usual three Mandelstam variables s, t and u, 
and is identified with the physical decay amplitude at the point s = M]^, t = u = M^. At 
the so-called soft-pion point (SPP), where the momentum of one of the two pions is sent 
to zero {s = u = M^, t = M]^), this amplitude is related to the K it amplitude, up 
to 0{M^) corrections. The problem is how to extrapolate the amplitude from the SPP to 
the physical point. The only working method proposed so far has been to use CHPT at 
tree level - using the one-loop relation does not solve the problem because a number of 
unknown low energy constants appear 0]. 

In the present letter we set up a dispersive framework for the K tttt amplitude. We 
show that by solving numerically the dispersion relations one can do the extrapolation in 
a controlled manner. The unitarity corrections due to rescattering of the pions in the final 
state, and those due to ttK (virtual) rescattering in the t and u channels, can be accurately 
accounted for by solving the dispersion relations. These effects, which also appear to one 
loop in CHPT, are not the only source of possible large corrections to tree level: two 
subtraction constants appear which may also suffer from large 0{M]^) corrections. The 
soft-pion theorem provides the means to determine one of the two subtraction constants, 
up to terms of order M^. The other constant (the derivative in s of the amplitude at the 
SPP) is unfortunately not yet determined with the same accuracy, and at present can be 
estimated only with tree-level CHPT. A better determination of this constant is the core 
of the problem. Once solved the K tttt amplitude can be obtained with substantially 
smaller uncertainties than at present. 

Truong 0, and more recently Pallante and Pich have stressed the importance of 
final state interactions in ^ vrvr, for the A J = 1/2 rule and e'/e, respectively. In 
estimating these effects they rely on the dispersion relation for the K ^ irn amplitude 
with the kaon off-shell. While the method provides a quick and simple estimate of the 
effect of final state interactions, it is not trivial to promote it to a systematic and rigorous 
calculation. The problems related to the formulation and the use of dispersion relations 
for an off-shell amplitude are discussed in a separate note ^ (see also 0). 
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2. We consider the amplitude^ 

(7r(pi)7r(p2)(/ = Q)\74'mK{q,)) =: T+(s,t,n) 
described in terms of the Mandelstam variables: 



s = (p^ +p2)^ t = [qi-pi] 



u 



{(ii-P2y 



(1) 



(2) 



related hj s + t + u = 2M^ + M]^ + g|, where q2 is the momentum carried by the weak 
Hamiltonian. From now on we set = (but ^2 7^ general). The physical decay 



amplitude is obtained by setting q2 = {s = M]^, t 



u 



Ml). 



Since the weak Hamiltonian has the quantum numbers of the kaon, and the pions are 
in an isospin zero state, the amplitude (|l|) is analogous to the t ^ u even combination of 
the KK —>■ nir scattering amplitude (the notation is borrowed from Ref. 0). Like in that 
case, one can show that if one neglects the imaginary parts of D waves and higher in all 
channels, then the analytic structure of the amplitude simplifies and it can be decomposed 
into a combination of functions of a single variable (for the Ktt scattering case see ||10|| ): 



[s,t,u) 



+ 



Mo{s) + ^ [Noit) + Noiu)] + ^ [Roit) + Ro{u)] 

s-u- ] N^{t) +is-t- 1 N,iu) 



t 



u 



(3) 



where A = Mj^ — Ml. Notice that the terms proportional to A^i drop out from the physical 
decay amplitude: 



A{K 



Mo(M|) 



+ 



Nq{MI) + 2R,{MI] 



(4) 



Each of the single variable functions appearing in Eq. (|^) is analytic in the complex 
plane except for a cut starting at 4M^ for Mq and at (M^- + M^^Y the remaining ones. 
These functions are defined to have the discontinuity on the positive real axis identical to 
that of a specific partial wave: Mq to the / = S-wave in the s channel, whereas in the 
t channel, Nq and Ni to the I = 1/2 5*- and P- wave respectively, and Rq to the 7 = 3/2 
S'-wave0. Below the inelastic threshold, the elastic unitarity condition for these functions 
reads 



discMo 




= sin5°(s 
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-< [Mo{s) 


+ Mo{s) 
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= sin^o''^! 


s] 


.,3/2 


Ro 


[s)+Ro{s) 



(5) 



^We discuss here only the AI —1/2 ampUtude - the A/ = 3/2 can be treated similarly. 
^We disregard the imaginary part of the 7 = 3/2 P-wave in the t channel because it is phenomenologi- 
cally very small and vanishes in the chiral expansion up to order 
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where Sq is the vrvr phase shift, whereas those with half-integer isospin are the ttK phase 
shifts. 

The hat functions denote contributions from the other channels coming in via angular 
averages (to be specified below), and are defined as 



Mo(s) 



{N,) + 2s\p\\q\{zN^) 



+ A\p\'\q\'{z'N,) + '^{{No)+2{Ro)) 



No{t) = {Mo)s + y{t){M^).,-r{t){zM^), + ^{4{Ro)^-{No)u) 

- ^-w{t){N,U-^v{t){zN,)^ + K\t){z'N,)^ , 

^iW = -^l{zMo)s + yit){zMi),-r{t){z^M,), + lmzRo)u- {zNo)u] 
r{t) I 3 

- lw{t){zN,)^ - ]v{t){z'N,)^ + lr\t){z'm)u 
o 4 o 

Roit) = {Mo)s-^yit){M^)s + ^r{t){zM^)s + ^i{Ro)u + 2{No)u) , 

+ ^w{t){N,U + ^v{t){zN^)^-^r\t){z'N,)^ , (6) 



where 



/Vf2 A 

iVi(t)=iVi(t)/t , S = Mi + M2, Si = S + M2, y(t) = Si-3t-^— , 
pMt) = V(l - (MK + M^y/t) (1 - {MK-M^)yt) , r(t) = (t - M2)p,;,(t) , 

= U - ^ , = - 4tSi + 5M^4 + S2 - . (7) 



The brackets () indicate angular averages defined as 

(z"X)(s) = ^ 1^'^ cizz"X(Si/2 - s/2 + 2|p| |q|z) 

{z''X)^{t) = - dzz''X{v{t,z)) v = s,u , (8) 
2 j-i 

where 

\P\' = 1-Ml \c\'='-(^l-MiJ. s + u = E,-t, s-u = ^ + r{t)z. (9) 

In the definition of the hat functions the function Mi appears. This is analogous to Mq, 
in the case of the 1=1 P-wave in the s channel, and is necessary to describe the process 
(ffl) in full generality, for all channels (including the t ^ u odd, 1 = 1 s-channel). It does 



not contribute directly to the physical decay process: its indirect (and small) contribution 
via the angular average in the dispersion relation is a consequence of crossing symmetry. 



3. If one is only interested in the low-energy region, neglecting the inelastic channels is a 
good approximation: then the solution of the dispersion relation for each of the functions is 
well approximated by the Omnes function times a polynomial |Tl| . It is therefore convenient 
to write the dispersion relation for the functions divided by the corresponding Omnes 
function (see for a detailed discussion of this point, although in a different framework), 
in the following form: 
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Niis) 
Rois) 
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n|(s) is the Omnes function [jlT|, defined as 
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(10) 
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(11) 



All functions are subtracted at s = with the only exception of Mq, where the sub- 
traction point So is left unspecified. In the following we use Sq = M^. The fact that 
only Mo depends on subtraction constants does not have any deep reason: the splitting of 
polynomial terms of T"*" between the various functions Mq, A^o.i and Ro is arbitrary, and 
we have merely used this freedom to remove them from the latter three. The final result 



does not depend on this choice |T2[ . All the dispersive integrals above have been cut off at 
energies Ai 2 and Ai 2 - numerical values will be given below. 



4. If the 7[7c and Kn phase shifts, the cutoffs Ai^25 and the subtraction constants a and b 
are given, the dispersion relations (0) can be solved numerically. Such a solution gives the 
amplitude T~^{s,t,u) at any point (provided it is far enough from the inelastic thresholds) 
of the Mandelstam plane, in particular at the physical point. The crucial new inputs here 
are the two subtraction constants: the phase shifts are known with sufficient accuracy, 
whereas the choice of the cutoffs is dictated by the inelastic thresholds. Before proceeding, 
we have to discuss how these two subtraction constants can be determined. If they could 
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be calculated with better accuracy than the physical amplitude itself, then this would 
represent for our method a clear advantage. 

For one of the two subtraction constants this is the case. The soft-pion theorem relates 
the amplitude at the SPP to the K ^ amplitude up to terms of order M^. We can 
therefore write 

" W/^^ - vr) = a + 1 [iVo(Mi) + 2Ro{MI)\ + 0{MI) , (12) 

which shows that a is indeed directly related to a quantity that is calculable (more easily 
than the decay amplitude itself), e.g. on the lattice. The relation (O) illustrates the 



strength of the soft-pion theorem: although the process involves a kaon, the relation is 
based on the use of the SU{2) symmetry, and therefore suffers from corrections of order 
only. 

The key of the problem is how to calculate b. This constant is related to the derivative 
in s of the amplitude at the SPP. The calculation of b requires the evaluation of the 
physical amplitude T'^ at an unphysical point, via analytic continuation. While this is easy 
to do with an analytical method like CHPT, it is practically impossible with a numerical 
method, like the lattice. However there is a Ward identity that relates this derivative to a 
Green function that is directly calculable: 

^T+(s, S - M2)|,=,,. = ^C{Ml Ml, Ml) + 0{Ml) , (13) 

where C{s,t,u) is an amplitude defined as: 



jdxe'P'%n{p2)\T'H\l,\Q)A^^{x)\K{q^)) = zp^B + iq>tC + iql^D , (14) 



where A^[x) is the axial current that couples to the pion removed from the outgoing state. 
By making the p2 momentum soft one can also derive a soft-pion theorem which relates 
the four-point function in Eq. (|1^) to a three-point function. Unfortunately the function 



C cannot be singled out from this relation. 

We are not aware of any attempts to calculate b. In order to illustrate our method we 
proceed by fixing 6 at a certain value and then varying it within a fairly wide range. To fix 
the central value and the range we use CHPT as a guide. At leading order, CHPT dictates 
the following relation between a and b: 

The size of the correction is at the moment unknown, but nothing protects it from being 
of order M|-: X = M|'/(167r^F|-)x, with x expected to be of order one. An explicit 



calculation in CHPT yields ||13|| : 



^ ^_^lnl + llni^-^(iV5 + 2iV7-Arg-4iV,o-4iVn) 
9 54 3 3 M2 3^ ^ ^ ^ ° ' 



^We have dropped the contribution coining from the weak mass term - more on this below. 
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E(GeV) 

Figure 1: The function \T~^{s, t, u) \ plotted vs. E = along the line of constant u = M^: 
the result of our numerical study for different values of X are compared to tree level CHPT. 



+ 2(iVi9-iV2o)-^(2iV2l+iV22 + 2iV23) + o(^^j , (16) 

where iVj = 167r^N[{Mx) are the renormalized low-energy constants introduced in |I4[ . 
Since we lack information on many of the constants, the CHPT calculation ([16|) does not 
allow us to do more than an order of magnitude estimate for b. 

5. In our numerical study we have used X = ±30%. The normalization of the amplitude 
is irrelevant here, and we have fixed it to T~^{M^, Mj^, M^) = 1. The vrvr phase shifts are 



taken from [|15|, with the scattering lengths determined in |T^, and the nK phase shifts 
from [|T^. For the cutoffs we have used Ai = 1 GeV, A2 = 1.3 GeV, and Aj = l.OSAj. Our 
results are shown in Fig. 0, where we have plotted |T+(s, S — s, M^)| versus s, comparing 
our numerical solution of the dispersion relations to the CHPT leading order formula. The 
latter is what has been used so far whenever a number for the K tctt matrix element 
extracted from the lattice has been given. Our treatment shows that large corrections 
with respect to leading-order CHPT are to be expected. One source of large corrections 
is the Omnes factor due to vrvr rescattering in the final state P, |]. The other potentially 
dangerous source is represented by X, the next-to-leading order correction to the relation 
(p!5|) between a and b. The latter could (depending on the sign) in principle double, or to 
a large extent reabsorb the correction due to final state interactions. The dependence on 
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X is well described by the following linear formula: 



l-^i^-^-^l 1.5(1 + 0.76X) , (17) 



I^Lo CHPT^;;^ _^ tttt) 



after having normalized both amplitudes to T~^{M^, M|^, M^) = 1. The evaluation of the 
uncertainties to be attached to the numbers in Eq. is in progress [0. At the moment, 



however, the main source of uncertainty is the fact that X is largely unknown. 

One of the outcomes of the present analysis is that the effects embodied in the functions 
A'o,! and Rq have turned out to be very small: if we drop them altogether, the numbers 
in Eq. (|l^) change from 1.5 to 1.4 and from 0.76 to 0.75. Notice that these effects are in 
principle of order M]^, as can be seen in Eq. (|T2p, and that they are not a priori negligible. 
On the other hand this result is very much welcome, because the size of these functions 
depends both on the irK phases (which are less well known than the tth ones) and on the 
choice of the cutoff A2, which may induce large uncertainties. 



6. The authors of have given a recipe to remove by hand the contribution of the weak 
mass term to the K ^ n transition. This step is necessary since at tree level the weak 
mass term does not contribute to the K ^ tttt amplitude. The procedure suggested in |^ 
involved the calculation of the K — * vacuum transition. In our framework this procedure 
is unnecessary. To show this it is useful to work with the tree-level CHPT amplitude (weak 



mass contribution only, see [0 for the notation) 



-ic^A 



Ml 



+ 1 



;i8) 



where we have restored ^ io show explicitly the presence of the kaon pole. At the 
physical point this amplitude vanishes. On the other hand, if we calculate this amplitude 
at the SPP, from the corresponding K ^ tt amplitude, we cannot get the kaon-pole term, 
and therefore would obtain a nonvanishing contribution at the physical point. Hence the 
cure proposed in [0]. The fact that the K ^ n amplitude does not contain the pole term 
does not contradict the soft-pion theorem: at the SPP, for = 0, the pole term is of order 
M^, which is beyond the accuracy the theorem guarantees. 

The framework proposed here relies on the presence of two subtraction constants, that 
we have related to the amplitude and its derivative at the SPP. If one determines both 
constants at tree level in CHPT (18), one finds that the subtraction polynomial vanishes 
at the physical point, without further ado. This cancellation does not happen if the 
subtraction constants are given to one loop - but it does not need to: as is well known [Q, 
the K TTTT amplitude does contain contributions proportional to C5 at one loop. 



7. In the present paper we have set up a dispersive framework for the K —>■ tttt amplitude 
that allows one to evolve the amplitude from the soft-pion point (where it is given by the 
K IT amplitude) to the physical point, taking into account all the main physical effects. 
As we have pointed out, this evolution is on safe ground only if a second input is made 
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available: the derivative of the amplitude at the soft-pion point which, to the best of our 
knowledge, has not been calculated so far. We have calculated this second subtraction 
constant to next-to-leading order in CHPT. Given the presence of unknown low-energy 
constants, we cannot use this expression for more than an order of magnitude estimate. 
Our numerical work, however, shows that the amplitude at the physical point depends 
strongly on the value of the slope at the SPP, see Fig. 1 and related discussion in the text. 
A nonperturbative calculation of the second subtraction constant b is necessary in order 
to obtain an accurate result with this method. We have provided a Ward identity which 
might be useful in this respect. 

Lattice calculations of the K tttt amplitude made so far rely on tree-level CHPT to 
relate the calculated K ^ ir matrix elements to the physical decay amplitude. The method 
proposed here improves this scheme by combining input from the lattice with dispersive 
techniques, thereby providing a fully consistent treatment of final state interactions in 
K ^ nil. Given the two subtraction constants, the dispersion relations can be solved 
numerically to good accuracy. Recently, a direct calculation of the K ^ nir matrix element 
on the lattice has been proposed in Ref. [|T^] - this method does not rely on CHPT. Other 
lattice methods, which also do not rely on the evaluation of the K ^ n amplitude had 



also been proposed previously Each of these methods presents different technical 

problems in its practical implementation and it is difficult to predict which one will 
lead to a reliable calculation of the K tttt amplitude. We hope that the present work 
will stimulate further efforts to calculate the subtraction constants a and b, either on the 
lattice, or by other nonperturbative methods. 
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